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Abstract 

We present a study of 49 single field, slow roll inflationary potentials 
in which we assess the likelyhood of these models fitting the spectral para¬ 
meters of the CMB radiation, namely the spectral index n s , its running a a , 
the running of the running (3 S and the ratio of tensor to scalar perturbations 
r, according to the currently most accurate determination of these spectral 
parameters given by the PLANCK collaboration. A double, partly redundant 
approach is followed: for most models we derive analytically bounds for the 
values of the spectral parameters that they can support, and for all of the mod¬ 
els we check numerically with a MATLAB program the spectral parameters 
that each model can yield for a very broad, comprehensive list of possible 
parameter and field values. 

The comparison of spectral parameter values supported by the models 
with their successive determinations by the PLANCK collaboration leads to 
contradictory conclusions: PLANCK2013+WP+BAO:ACDM+/+a s disfa¬ 
vours all of the models with confidence at least 93%, conversely the data 
provided by PLANCK2015+TT+lowP without running of the running allows 
back most of the models, but taking into account the running of the running 
again disfavours 39 of the 49 models with confidence at least 92.8%, and 5 
more because of the amount of expansion e-folds supported. 

We identify a bias in the method of determination of the spectral paramet¬ 
ers currently used to reconstruct the power spectrum of scalar perturbations 
that can explain these contradictory conclusions. The solution to this prob¬ 
lem is likely to determine the fate of the inflaton. 


1 Introduction 

The PLANCK2013 data |[Q showed that the spectral index of scalar perturbations 
has an expected negative running, whose modulus was in general one order of mag¬ 
nitude greater than the theoretical values provided by slow roll inflation, which 
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could be used to test some theoretical slow roll inflationary models. The aim of 
the present work is to present an analytic and numerical study, of the 49 theoretical 
models of single field, slow roll inflation that are described in 0 (see its Table 
1), that shows that the constrains of the spectral index, its running and the ratio of 
tensor to scalar perturbations provided by PLANCK2013 combined with several 
data, disfavours all of these models, severely in most cases. The most familiar 
ones such as Hill-top, Natural, Plateau and Monomial potentials could be disreg¬ 
arded according to these PLANCK2013 data due to the small value, in modulus, 
of the theoretical value of the running provided by slow roll inflation. In fact, from 
PLANCK2013+WP+BAO: ACDM+r+o iS data we show that the deviation from the 
theoretical value of the running, namely D, obtained from the majority of that mod¬ 
els to its expected observational value is larger than 1.6cr, and in modulus all the 
theoretical values are smaller than the modulus of the expected observational one, 
which means that these models lies outside of the 94.5% C.L. Further numerical 
study has shown that for all of the models either a s or n s or r lie outside of the 
93% C.L. 

Our study, initially performed from PLANCK2013 data, could suggest that 
single scalar inflationary theories must be replaced by multiple fields theories, by 
other ones with a breakdown of the slow roll phase 0 or by reconstruction tech¬ 
niques 0- Another completely different proposal is to abandon the inflationary 
paradigm in favour of the Matter Bounce Scenario (5[ |6l [71 where the Big Bang 
singularity is replaced by a non-singular Big Bounce, which at this moment, con¬ 
stitutes a promising alternative to the slow roll inflationary paradigm, without the 
characteristic inflationary flaws, such as the initial singularity which or the fine- 
tuning of the degree of flatness required for the potential in order to achieve suc¬ 
cessful inflation li8l . 

Fortunately for single field slow roll inflation, the new PLANCK2015 obser¬ 
vational data 0, reduce the modulus of the running one order, which allows the 
viability of the majority of single scalar field slow roll inflationary models. For ex¬ 
ample, using Planck2015+TT+lowP:ACDM+? +Q! s , where n s = 0.9667 ± 0.0066 
and a s = —0.0126^QQQg® at lcr C.L. (see table 4 of 0), and the conservative con¬ 
strain r < 0.25 (see figure 6 of 0), for the majority of tested potentials one has 
D = 1.1a, which means that these models are only disfavored at 86% C.L. or less. 
Moreover, if one introduces lensing, then D will be reduced to be lower than la, 
and thus allowing single scalar field slow roll inflation. 

However, dealing with the running of the running, namely j3 s , when one con¬ 
siders the PLANCK2015 TT+lowP (resp. PLANCK2015 TT,TE,EE+lowP) model 
/3 S = 0.029^0 (resp. f3 s = 0.025 ± 0.0013) (see (19) of 0), the same meth¬ 
odology that we have applied in the previous cases show that the measured value 
of f5 s is incompatible with the theoretical value provided by single field slow roll, 
because its measured value is too large for a magnitude that depends on third order 
on the slow roll parameters, disfavouring almost all slow roll inflationary models 
studied, 39 out of 49 at 92.8 % C.L. or more. 
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The main lesson of this work, is the importance of the observational meas¬ 
ures provided by PLANCK or other teams, in order to check the viability of single 
field slow roll inflation. Only after a precise determination of their expected ob¬ 
servational value and the corresponding deviation we will able to determine which 
inflationary models could depict correctly our Universe. 

The problem of the reliability of the observational measures of spectral para¬ 
meters of the CMB starts with the above paradox of the the oscillating conclusions, 
from ruling out every model to allowing most of them and back, that the successive 
PLANCK2013 and PLANCK2015 determinations of spectral parameters support. 
We find that these oscillating conclusions can be explained by a bias in the method 
of fitting spectral values of the power spectrum to the observations. We believe that 
new data, and the addressing of this bias, will reduce significantly the observable 
value of the running of the running and may possibly allow back the viability of 
single slow roll inflation, in the same way as has happened with the early values 
of the running obtained by PLANCK2013 and its drastic decrease according to 
PLANCK2015. 

The units used in the paper are: h = c = 87 tG = 1. 


2 Slow-roll parameters 


In slow roll inflation (see ifTOll for a review of inflation) the commonly used first 
order parameters are: 


e 
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At the first slow roll order, the spectral index of scalar perturbations and its 
running are given by 

n s — 1 = 2?7 — 6e and a s = 16 erj — 24e 2 — 2£, (2) 


where the second order slow roll parameter 
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has been introduced. 

Moreover, in inflationary cosmology, the tensor/scalar ratio, namely r, is re¬ 
lated with the slow roll parameter e, via the following consistency relation r = 16e. 

The other important parameter that we will use in this work is the running of 
the running j3 s = , given, in the slow roll approximation, by ifTlITOl 

p s = - 192e 3 + 192e 2 r ? - 32 er] 2 - 24e£ + 2r/£ + 2(, (4) 
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where we have introduced the third order slow roll parameter 
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3 Analytical fitting of the parameters 

3.1 PLANCK2013 data: the running 

In this section we will take into account the PLANCK2013 constrains on n s , a s 
and r which are model dependent (see table 5 of (TJ), and the fact that more than 
50 e-folds are needed to solve the horizon and flatness problems of GR. If one 
does not consider the running and makes the analysis in the plane (n s ,r), then 
PLANCK2013 data shrink the space of allowed standard inflationary models, pre¬ 
ferring potentials with a concave shape (V w < 0) IDQ]. But, it is the combination 
of the three data ( n s ,a s ,r ) what rules out all the standard slow roll inflationary 
models 0. 

Effectively, for instance, we will consider in detail the ACDM+r + a s model 
from PLANCK2013 combined with WP and BAO data, which gives the following 
results (see table 5 of [I]).: 

n s = 0.9607±0.0063, r < 0.25 at 95% C.L. and a s = — 0.021±g;g^g. 

In slow roll inflation, a simple calculation leads to the relation 

= 7 , ('As - l)r + i’ 2 - 2£. (6) 

And thus, assuming that the potentials we consider have spectral index n s at 
less than 2 a deviations and satisfy the conservative bound r < 0.32 (see figure 

4 of HI). the minimum of the function n s — 1 )r + ^r 2 reached at the point 
(n s = 0.9481, r = 0.1384) is greater than —0.0018, what provides the bound 

> -0.0018 - 2£, (7) 

meaning that plateau potentials such as V(ip) = Vo ^1 — 'j (Hill-Top Inflation 
(HTP)) El, V(<p) = Vo (l - ^) 2 with |<p| < /j (Double-Well Inflation (DWI)) 

lfl3l or V(ip) = Vo ^1 + cos (^j^j (Natural Inflation (NI)) lfl4ll . when one con¬ 
siders values of the running at lcr C.L., are disfavoured by PLANCK data because 
for all of them £ < 0. In fact, the deviation from the theoretical value of the running 
to its expected observable value, namely 2), is larger than 1.6<r. 

A distance larger than 1.6<r is also obtained for the potential that leads to Ex¬ 
ponential SUSY Inflation (ESI) U(<^) = Vq (1 — e ~ w ) lfl5l and for Power Law 
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Inflation (PLI) whose potential is given by V{p) = Vbe PV> lfT6ll . because in this 
case one has £ = p i 2 * , that is, 


16 


3 r 2 


= - (n s - 1) + — r - ~{n s - 1) = —— > -0.0018, 


32 


( 8 ) 


where we have evaluated a s , as a function of n s and r at the absolute minimum, 
namely n s = 1 — ^ with r = 0.1384, in the rectangle [0.9481, 0.9733] x [0,0.32]. 
To be more precise, since power law inflation has no running, the distance to the 
main observational value is 1.7cr, because in this case one has the bound a s > 
-0.00045. 

Dealing with the Kaller Moduli Inflation I (KMII) iTTll given by the potential 
V{p) = Vo (1 — ape~ v ), with Higgs Inflation (HI) |fT8l where the potential is 

V{p) = Vo ^1 — e and for Open String Tachyonic Inflation (OSTI) lfl9ll 

with potential V {ip) = —Vq In , one always has £ < rj 2 and one obtains 
the same conclusions as in ESI, i.e., 2) > 1.6c. 

When one considers Large Field Inflation (LFI) given by the monomial poten¬ 
tial V{ip) = Vq p> v with p > 1 Il20l . one obtains 


n, - 1 = - 


Pip + 2 ) 

<P 2 


and 
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which means that p must be positive in order to have an spectral index with a red 
tilt and a negative running. As a first consequence, Inverse Monomial Inflation 
(IMI) Ell is disfavored. Second, from the relation 


a, =- n s - l) 2 > -\{n s - l) 2 > -0.0018, (10) 

we conclude that 2) > 1.6cr. And for Radiation Gauge Inflation (RGI) Il22l l. whose 
potential is given by V(ip) = F 0 ^2 one has £ = V 2 < because 

i ncreases as a function of ip. Then, evaluating at n s = 0.9481 one has 


ot s = —^-(n s — l) 2 > —0.0018, (11) 

O 

giving as a result 2) > 1.6 <j. 

For potentials such as: V(<p) = Vo ^1 — ^ (Brane Inflation (BI)) lf23l . 

and V{ip) = Vo ^1 + ^ (Dynamical Supersymmetric Inflation (DSI)) f24l . 

since one has £ = p 2 , one can obtain the following exact formula 
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The minimum of a s is obtained at n s = 1 — ^ with r = 0.1384, and thus, 
inserting this expression in ([T2]), one gets 


3 r 2 

> -0.0018, 


(13) 


which is incompatible with the running provided by PLANCK at la C.L., because 
D > 1.6cr. 


For general hill-top potentials such as: V(p) = V'o ( 1 — 


(Small Field 


Inflation (SFI)) lfl2l and V(<p) = Vb (1 + ( ^) ) (Valley Hybrid Inflation (VHI)) 

Pz 
p- 


with p > 3, since one has £ = one can obtain the following exact 


formula 


CX.Q — 


(p + 2 )r 


^ (p-2)(n s -l) 2 

— 1) + iej- 2(p — i) ' 


(14) 


8 (p ~ 1) 

Since the minimum of a s is obtained at n s = 1 — with r = 0.1384, one also 


has 


3r 2 

a s > -> -0.0018, 

32 “ 


(15) 


giving D > 1.6(7. 

However when one deals with Arctan Inflation (AI) Il26ll with potential V(cp) = 
Vo ^1 — ^arctan where one has £ < ij// 2 . The absolute minimum is 

reached at the point (n s = 0.9481, r = 0.32), leading to the constrain 


(n s — 1 )r 3 r 2 3 (n s — l) 2 

a s = - - -—-> - 0 . 0022 , 

16 256 4 “ 


(16) 


which means that the deviation to the expected observational value is greater than 
1.56(7. 

In the case of Loop Inflation (LI) |[27l with potential V(ip) = Vo(l + a In ip), 
one has £ = 2 rj 2 , leading to the constrain 


a s = —^(n s - l)r - ^-r 2 - (n s - l) 2 > -0.0034, 


(17) 


when one evaluates at the point where a s reaches its absolute minimum, namely 
(n s = 0.9481, r = 0.32). Consequently, D > 1.46(7. 

For Mixed Large Field Inflation (MLFI) lf28l with potential V(ip) = V()^ 2 (l + atp 2 ), 


one has £ = 2 (i+ 2 a<p 2 ) r - l* 16 ot h er hand, the relation 


r = 16e = 


32 /1 + 2 
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leads to the bound 1+ ^ aip2 < qqqq < gfe, and thus, £ < Finally, evaluating 
at the absolute minimum (n s = 0.9481, r = 0.276) we can conclude 

1 3 3r 2 

a s > ~(n s — 1 )r 4-r 2 =-- > —0.0036, (19) 

2 v ' 64 64 “ ’ v 7 

and thus, D > 1.45<r. 

Finally, in Witten-O’Raifeartaigh Inflation (WRI) |29l with potential V(tp) = 
Vo In 2 one has £ < | rj 2 , leading to the constrain 

11 9 

a s > — — Rn s — l)r-r 2 - (n s — l) 2 > —0.0040, (20) 

32 v s ' 512 8 V ' “ 

when one evaluates at the point where a s reaches its absolute minimum, namely 
(n s = 0.9481, r = 0.32). Consequently, £) > 1.42<r. 

To end this section a remark is in order: First of all, it is important to real¬ 
ize that we have studied analytically only 25 of the 49 models provided by d- 
Secondly, our analytic results have been obtained bounding the minimum of a s in 
the rectangle 

R = {(ra a , r) : 0.9481 < n s < 0.9733, 0 < r < 0.32}. 

However, the running could be parametrized with only one independent variable, 
for instance, the scalar field ip, which means that the bounds obtained in the rect¬ 
angle could be improved, because one only needs to find a bound inside a curve 
inside the rectangle R. The problem to perform this calculation analytically is 
that this curve can only be obtained explicitly for a few potentials. The numerical 
calculations of Section[4]are needed to improve the analytical ones. 

3.2 Analytical results from PLANCK2013: A CDM+/ +o s model com¬ 
bined with other data 

The same kind of results could be obtained from other models with running. In 
fact, assuming that the potentials we choose have an spectral index n s at less than 
2(j deviations and also they satisfy the conservative constrain r < 0.32, we have 
summarized the results in four Tables: Table 1 contains the description of three 
models and Table 2 (WP+high-/:), Table 3 (WP) and Table 4 (WP+lensing) contain 
the deviation of the running, for the potentials we have analytically studied, for 
each one of the models. Again, the computed likelihoods reflect that the values of 
the running a s predicted by the potentials lie all in the same tail of the Gaussian 
distribution. 
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Determination 

r 

n s 

OL s 

PFANCK2013+WP+ high -l 

< 0.23 

0.9570 ± 0.0075 

—0 022 +0 ' 011 
o.uzz_ 0 010 

PFANCK2013 +WP 

< 0.25 

0.9583 ±0.0081 

-0.021 ±0.012 

PFANCK2013 +WP+lensing 

< 0.26 

0.9633 ± 0.0072 

-0.017 ±0.012 


Table 1: PLANCK2013 estimations of spectral parameters, without running of the 
running (ID). 


Potential 

Running deviation 

Disfavored at 

PFI 

> 2 a 

97.75 % C.F. or more 

HTI, DWI, NI, ESI, KMII, HI, PSNI 



SFI, FFI, VHI, DSI, BI, OSTI, RGI 

> 1.8(7 

96.4 % C.F. or more 

AI 

> 1.7a 

95.55 % C.F. or more 

FI, MFFI, WRI 

> 1.6a 

94.5 % C.F. or more 


Table 2: PLANCK2013+WP+high-(:ACDM+r+a. s model. 


3.3 PLANCK2015 data: the running almost vanishes 

The new observational data provided by reproduced in Table 5 reduces one or¬ 
der or more (depending on other data such as BAO, lensing and lowP), the modulus 
of the expected observational value of the running. 

This drastic reduction of the running in modulus is what allows the viability 
of the potentials disregarded from PLANCK2013 data. Effectively, using for in¬ 
stance, the PLANCK2015 TT+ 1 o w P+B AO: A CD M+r+a s data, and assuming that 
the potentials have spectral index n s at less than 2er deviations, and also satisfying 
the conservative bound r < 0.32. Then, for the simplest potentials, analytically 
one can show: 

1. For LFI, one gets the following bound 

|a s | < ^~p( ns ~ X ) 2 - ( n * - 1 ) 2 

which means T> < 1.2a, and thus, disfavoring 
86.45 % C.F.. 

2. For PFI, which has no-running, one obtains D < 1.4a, and thus, disfavoring 
the potential less than the 91.95 % C.F.. 

3. For HTP, ESI, BI, DSI, SFI and VHI the minimum and maximum of a s are 
obtained respectively at n s = 1 — ^ with r = 0.1101 and r = 0.0642. And 


< 0.0017, (21) 

the potential less than the 

















Potential 

Running deviation 

Disfavored at 

PLI 

> 1.7a 

95.55 % C.L. or more 

HTI, DWI, NI, ESI, KMII, HI, PSNI 



SFI, LFI, VHI, DSI, BI, OSTI, RGI, AI 

> 1.5a 

93.93 % C.L. or more 

LI 

> 1.4(7 

91.95 % C.L. or more 

MLFI, WRI 

> 1.3(7 

90.3 % C.L. or more 


Table 3: PLANCK2013+WP:ACDM+r+a,. model. 


Potential 

Running deviation 

Disfavored at 

PLI 

> 1.4(7 

91.95 % C.L. or more 

HTI, DWI, NI, ESI, KMII, HI, PSNI 



SFI, LFI, VHI, DSI, BI, OSTI, RGI, AI 

> 1.2(7 

88.5 % C.L. or more 

LI, MLFI, WRI 

> 1.1(7 

86.45 % C.L. or more 


Table 4: For PLANCK2013+WP+lcnsing: ACDM+r+a., model. 


thus, one gets —0.0012 < a s < —0.0003, what implies T) < 1.4(7, and thus, 
disfavoring the potential less than the 91.95 % C.L.. 

The rest of potentials in the examined list fit even better the values of the spec¬ 
tral parameters of Table 5 for some choice of parameter or field values. 

3.4 PLANCK2015 data: the running of the running 

The last PLANCK2015 data about the running a s and its running f3 s = are 
reproduced in Table 6. 

These values are similar to the ones of PLANCK2013, reproduced in Table 7. 
These results contradict single field slow roll inflation, because in that case, the 
running a s is second order in the slow roll parameters, and its running is given by 
Eqs. <(4]), Q, which make d s a third order parameter, while the values determined 
by PLANCK place /3 S in a higher order of magnitude than the running a s itself. 
Moreover, disregarding the running of the running, the running is negative while 
taking into account it, the running becomes positive. This seems a signature of 
the problem that suffers the method used to reconstruct the power spectrum of 
scalar perturbations from observational data: the value of the coefficients in the 
Taylor series of the power spectrum logarithm function could suffer a bias. We 
will address this question in Section 5. 
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Determination 

r 

n s 

OL s 

PLANCK2015 TT+ lowP 

< 0.180 

0.9667 ± 0.0066 

— 0 0126 +0 ' 0098 
u.uizo_ 0 0087 

PLANCK2015 TT+ lowP+lensing 

< 0.186 

0.9690 ± 0.0063 

— 0 0076 + °' 0092 
u.uu(o_ 0 0080 

PLANCK2015 TT+ lowP+BAO 

< 0.176 

0.9673 ± 0.0043 

-0.0125 ±0.0091 

PLANCK2015 TT, TE, EE+ lowP 

< 0.152 

0.9644 ± 0.0049 

-0.0085 ± 0.0076 


Table 5: PLANCK2015 estimations of spectral parameters, without running of the 
running (0). 


Determination 

n s 

Q. s 

Pa 

PLANCK2015 TT+lowP 

0.9569 ±0.0077 

0 011 + 0 - 014 
u ' uii -0.013 

0 099+ 0 - 015 
u.uzy_ Q 016 

PLANCK2015 TT,TE,EE+lowP 

0.9586 ± 0.0056 

0.009 ±0.010 

0.025 ±0.013 


Table 6: Determinations of the spectral parameter values, with running of the run¬ 
ning by PLANCK2015 ((19) of @) 


To show the improbability of the observed value of j3 s analytically for all the 
potentials that appear in 0 is very involved due to the increasing complexity of 
the formulas ([4]), ([5]) for the new parameter /3 S . This work thus follows a double 
approach, combining analytical calculations for LFI, LI (with a > 0), VHI, SFI, BI 
and DSI (with p and even number or p > 0), HTI, ESI and PLI, with a numerical 
analysis to which we subject all the potentials in 0. 

Let us start with the analytical approach. In the case of LFI one has 


n, - 1 = - 


p{p + 2 ) 


Ota - 


2 p 2 (p + 2 ) 


p- 




and (3 S = — 


8 p 3 (p ± 2 ) 


c p v 


( 22 ) 


and thus, 


P. = ~ 


(p± 2) : 


;(n s - l)' 


I A, I < gK - 1| 3 - 


(23) 


Then, for the PLANCK2015 TT,TE,EE+lowP data, considering n s at 2cr C.L., 
and thus, after inserting n s = 0.9474 in (23), one obtains the bound (3 S < 0.00013, 
which means that the deviation from the theoretical value of the running of the 
running to its expected observational value is larger than 1.9cr. A deviation larger 
than 1.9(7 is also obtained from PLANCK2015 TT+lowP data, meaning that LFI is 
completely disfavored by PLANCK2015 data. 

In the case of LI (with a > 0), VHI, SFI, BI and DSI (with p and even number 
or ip > 0), HTI, ESI and PLI, a simple calculation shows that 6 S < 0 which means 
that the deviation from the theoretical value of the running of the running to its 
expected observational value is also larger than 1.9 cj. 
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Determination 

n s 

cx s 

13s 

PLANCK2013+WP+BAO 

0.9568t°;°! 

o ooo + 0013 

U.UUU_Q 016 

0 01 7+o.ot6 
U - Ui ' -0.014 

PLANCK2013+WP+high-^ 

0.947618:811 

0 001+ 0013 
u.uui_ 0 014 

0 022 +0 ' 016 
u.uzz_ 0 013 

PLANCK2013+WP+lensing 

0.957318:8^ 

0 OOfi +a015 
u.uuo_ 0 014 

0 01 Q+ a018 
u.ury_ 0 014 

PLANCK2013+WP 

0.951418:8^ 

o ooi +0,016 

u.uui_ 0 014 

0 020 +ao16 
u.uzu_ 0 015 


Table 7: Determinations of the spectral parameter values, model ACDM+a s +/3 s , 
by PLANCK2013 (table 5 of Q) 

4 Numerical fitting of the parameters 

4.1 Methodology 

Let us describe the numerical tests that the authors have applied to all single field 
inflationary models from the list of (2J. These tests have been built into a MAT- 
LAB program that takes as input a list of potentials V(ip) and values for spectral 
parameters in the list r, n s ,a s , /3 S , and asseses the likelihood of each model in the 
list to fit the values of the spectral parameters using a 95% confidence limit for the 
value of r, and assuming Gaussian distribution for the values of n s , a s , (3 S . 

For each cosmological model, a broad range of possible values for the paramet¬ 
ers on which it depends has been determined following |0. A test list of values for 
each parameter has been selected, covering in a dense, equispaced fashion finite in¬ 
tervals of possible values for the parameter, and approaching with log-equispaced 
values every finite or infinite limit value for the parameter. 

After allowing simplifications induced by rescaling, if a model still depended 
on more than one parameter all possible combinations of values for each parameter 
were tested. Table [8] lists the selected parameter values for each model. 

Also separately for each model, the range of possible values for the inflaton 
field was determined taking into account whether the model admitted values of 
the field with any sign, or only positive values, and further peculiarities such as 
periodicity of the model and the applied rescalings. The considered range p/} 
for field values in each model is also listed in Table [8j 

The MATLAB software developed by the authors, for each model V (p) and 
choice of value of the parameter(s) on which it depends, takes an equispaced mesh 
of values in the range Pf] °f possible values of the field in this model. This 
mesh is taken increasingly fine, currently up to step A ip = 2 • 10~ 4 . 

The subintervals in the range of field values for which the potential satisfies 
V(p) > 0 are numerically determined over the selected mesh, and each interval of 
positive values of the potential for the selected values of the model parameters is 
considered as a case , which thus consists of: 

• a candidate theory with a given potential V(ip), 
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• a specific choice of parameter values for V, 

• and a range of values [<po,<pf\ of the inflaton field ip such that V > 0 on 
them. 

The numerical test for each case consists in meshing the interval of field values 
with a uniform step (of size A p = 2 • 10 4 for the results reported in this work), 
computing the spectral parameters r, n s ,a s , (3 S for each value of the field <p in the 
mesh using the formulas of Section [2] and symbolic derivation of the potential V 
to produce the derivatives V ^,..., V ww , and then applying successive filtering 
criteria to determine which values of the field ip fulfill simultaneously all of them, 
thus allowing the model in this particular case to fit the spectral measured data for 
which the model is tested. 

A determination of a set of values for the spectral parameters provides their ex¬ 
pected values < n s >, < a s >, < j3 s >, and their standard deviations a Us , a as , og s . 
In the case of spectral parameters without running of the running (3 S , this is re¬ 
placed by the tensor/scalar ratio r, which is conservatively estimated to have a 
value r < 0.32 (see lffll.f9l'). 

The applied filters in the case of spectral parameters without running of the 
running consist in looking for the values of the field ip such that: 

1. r < 0.32, 

2. | n s (ip)- < n s > | < 2 a ng , 

3. |a s (</?)- < a s > | < 1.6<7 as , 

A model not passing the first filter is ruled out with 95% C.L., a model not passing 
the second filter is ruled out with 95.5% C.L., and a model not passing the third 
filter is ruled out with 94.5% C.L., because the values of a s provided by the models 
are located uniformly in the same tail of the Gaussian distribution (namely, from 
the negative expected value towards zero) for the spectral valuations to which we 
have applied our test. 

The applied filters in the case of spectral parameters with running of the run¬ 
ning consist in looking for the values of the field ip such that: 

1. e(ip) < 1, 

2. |n s (<^)— < n s > | < 2a na , 

3. |a s (<£>)- < a s > | < 2a as , 

4. \/3 s (ip)~ < (3 S > | < 1.8cra s , 

A model not passing the first filter (equivalent to asking for r < 16) is ruled out 
with 95% C.L., a model not passing the second or third filter is ruled out with 
95.5% C.L., and a model not passing the fourth filter is ruled out with 92.8% C.L. 
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In either situation with or without running of the running, a case (choice of 
model, parameter values, and interval of values for the field) is considered possible 
only if there exists some value of the field c p in its range such that it satisfies sim¬ 
ultaneously all of the filtering conditions. The filters are not independent, thus if a 
case does not satisfy all of the filtering criteria simultaneously for any field value ip 
in its range, it is disproved with the confidence level of the strongest filter it fails. 
A model such that it is disproved in any case (i.e., for any choice of parameters 
and range of field values) is regarded as disproved with the lowest confidence level 
with which any of its cases is disproved. 

Conversely, if a case has values of the field ip satisfying all of the filtering cri¬ 
teria, the minimal values among them of the distances n s (p)— < n s >,a s (ip ) — < 
a s >, l3 s (ip)— < 13 s >, expressed in standard deviations, are logged. The max¬ 
imum among these distances provides the confidence level to which the particular 
case has been disproved. In the case of computations without running of the run¬ 
ning such that all of the values of a s provided by the models are found consistently 
in the same tail of the Gaussian distribution then - likelyhood is assessed taking into 
account only this tail. 

In these cases satisfying all filters, i.e. that are not disproved with C.L. at least 
92.8% the testing software also looks for values of the field tp e such that e(p e ) — 1 , 
and using them as endpoints of the inflationary phase, computes the number of e- 
folds of inflation for any choice of tp in the case, by integrating numerically with a 
trapezoidal rule 



(24) 


The result N(ip) ranges over the number of e-folds of expansion that the case 
supports for the field values p satisfying all of the filtering criteria. Its minimal and 
maximal values are logged, as they are the limit values for the number of e-folds 
of expansion that the case can support. 


4.2 Numerical results 

Single-field inflaton models were exhaustively studied in ||2l, from which we take 
the list of models and parameters to be numerically tested. Table [8] adapted from 
Table 1 of (2|, presents each model’s potential, the range of values of the paramet¬ 
ers for which it has been tested, and the range of values of the inflaton field over 
which it has been tested. 


Table 8: Models from [2] numerically tested. 


Name V ((p) Parameter values Field values 


HI 

Vo (l - e“V / 573 V’Y 


[-40,40] 

RCHI 

Vo(l-2 e -V2/3^ + _|2_^) 

Aj : [linspace(-100,100,120),linspace(-3,3,200)] 

[-10,20] 

LFI | V 0 (<p>) p | p: linspace(0.5,20,60) 

MLFI 

Vq<p 2 1 + onp 2 

a .: [linspace(-10,100,61 ),linspace(-0.1,0.1,120)] 
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RCMI 

Vo M 2 

1 — 2 o:y ) 2 In (c/ 2 ) 


a : [Iinspace( 1 e-4,1.5,30), 10.Alinspace(-14,-5,20)] 


RCQI 

Vo M* [1 - a In O)] 

o : [linspace( 1 e-2,10,120), 10. A linspace(-6,-2.5,10)] 


NI 

Vo 

1 + COS (y ) 


/: 1 

[ 10 _ 5 ,x] 

ESI 

Vo (l-e-^) 

q : [linspace(0.1,10,120),linspace( 1 e-5,0.099,60)] 


PLI 

Voe-*** 

o : [linspace(0.1,10,120), 10. Alinspace(-6,-1.5,15)] 

[-40,40] | 

KMII 

Vq ^1 — aipe~ 

o : [linspace( 1 e-2,10,60), 10. Alinspace(-6,-2.5,10)] 


HF11 

V 0 (1 + A 1V ) 2 1 2 3 ( 1+ % v ) 


A\ \ linspace(le-3,40,180) 

[-40,40] 

CWI 

Vo 

1 + “(§) ln (§) 


Q : [ 10. A linspace(-6,-2.5,12),linspace( 1 e-2,10,60)] 


LI 

Vo [1 + a In (<^)j 

or. [linspace(-0.3,0.3,60),-10.Alinspace(-1,0,10), 10.Alinspace(-1,0,10)] 


Rpl 

Voe -2y/m* _ 1 | 2 *>/C2*»-1> 

p\ linspace(0.25,10,60) 


DWI 

V 0 

(—) 2 - 1 
\ V’O ) 

2 

vo- 1 

[10 —4 , SO] 

MHI 

Vo 

i-sechf^) 


p: 10 

0 

0 

1 

0 

RGI 


ol : [linspace( 1 e-1,10,30), 10. Alinspace(- 6 ,-1.5,12)] 


MSSMI 

Vo 

fje.'j 2 -2 (j^ 6 + i 

V V’o / $ \<P 0 / 5 \ / 

J 

<p 0 \ [le-7,le-3,l] 

[10 —4 ,80] 

RIPI 

Vo 

M 2 _4M 3 ,iM 4 

V ‘PO / 3 V <P0 J 2 V ‘po / 

J 

(Pq: [le-7,le-3,l] 

[10 —4 , SO] 

AI 

Vo 

1 — y arctan 


p: [le- 2 , 1 ] 

[-40,40] 

CNAI 

Vo 

3 — (3 + ck 2 ^ tanh 2 ( -y= 


O ’: [linspace( 1 e-2,20,120), 10. A linspace(-6,-2.5,10)] 


CNBI 

Vo 

(3-o 2 )tan 2 (^¥>) -3 


o : [linspace( 1 e-3,5,40), 10.Alinspace(-7,-3.5,12)] 


OSTI 

-Vo(^) 2 ' n 

( —) 2 
V^o/ 


vo- i 

[ 10 “ 6 , 1 ] 

WRI 


vo- 1 

[ 10 — 4 , 10 ] 

SFI 

M 1 -(£)'] 

p: 1 

p : linspace(0.5,10,20) 

[ 10 - 4 , 1 ] 

II 

v 0 (v - Vo)~ P - Vo ^ (v - vo)~^~ 2 

VO - 0 

/ 3 : [linspace(0.1,10,3 l),linspace(20,50,3)] 


KMIII 

Vo 

4 

1 — a:(c^) 3 exp 



a: 10.Alinspace(-3,12,46) 

/3: 10.Alinspace(-3,12,46) 

[ 10 — 4 , 10 ] 

LMI 

Vo (¥>)“ exp [—PM'*] 

j3 : [linspace(0.1,20,40), 10. Alinspace(-4,-1.5,10)] 

7 : [linspace(le-3,2,30),10.Alinspace(0.5,2,4),10.Alinspace(-6,-4,3)] 


TWI 



¥>0 : 1 

A: linspace(0.001,8,120) 


GMSSMI 

Vo 

f— ) 2 - i a (— l 6 +t (— ) 10 

]_ 

<£>0 : 10.Alinspace(-2,0,3) 

0 : [linspace(le-2,2.5,120),10.Alinspace(0.5,1.5,3),10.Alinspace(-4,-2.5,6)] 


GRIPI 

Vo 

f-t) 2 - f (A) 4 

V ‘PQ J 3 \ Vo J 2 \ <PQ ) 


ipo'. 10.Alinspace(-2,0,3) 

a : linspace( 0 . 1 , 10 , 120 ) 

[ 10 — 4 , 20 ] 

BSUSYBI 

V 0 

7 : linspace(le-5,2,200) 

[-40,40] 

TI 

Vo (l + cos ^ + a sin 2 £) 

p: 1 

0 : [linspace(0.01,3,80), 10.Alinspace(-4,-2.5,6), 10.Alinspace( 1,2,4)J 


BEI 

V 0 expj(-Aip) 

/3 : [linspace(-5,5,60), 10. A linspace( 1,2,3), 10. A linspace(-4,-2,3), 

-10. Alinspace( 1,2,3),-10. Alinspace(-4,-2,3)] 

A: 1 

[- 100 , 100 ] 

PSNI 

Vq [l + a In (cos ^ J 

0 : [linspace(0.1,10,60), 10. Alinspace(-4,-1.5,12)] 

f -1 

1 

0 

1 

(M 

1 

O 
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NCKI 

Vo [l + Q In (ip) + p (ip) 2 ] 

or. 10.Alinspace(-7,0,16) 

linspace(-l 0,10,80) 


CSI 


a: [linspace(0.1,5,100), 10.Alinspace(-4,-1.5,8), 10.Alinspace( 1,2,4)] 

[-40,40] 

01 


vo'- 1 

cx: [ 10.Alinspace(-7,-2,18),linspace(0.03,1,40)] 


CNCI 

Vo (3 + a 2 ) coth 2 (^v) - 3 

a: [linspace(0.1,5,40), 10. Alinspace(-7,-1.5,12), 10. Alinspace( 1,3,5)] 


SBI 

Vo {l + [-a + /31n(<«o)] (<p) 4 } 

ck: 10.Alinspace(-8,0,27) 

(3: 10.Alinspace(-8,0,27) 


SSB1 

V 0 [l + a (v) 2 + p M 4 ] 

ol : [-10.Alinspace(-5,2,24), 10.Alinspace(-5,2,24)] 

(3: [-10.Alinspace(-5,2,24), 10.Alinspace(-5,2,24)] 

o 

1 

to 

o 

IMI 

V 0 (v)~ p 

p: linspace(0.5,10,40) 


BI 

V 0 

1 

1 

•e_ 


p: [linspace( 1,10,37), 10. Alinspace(-1 ,-0.33,3)] 

p: [le-4,1] 


RMI 

V 0 [l-f (-i+ln^)* 2 ] 

vo - 1 

c: [-linspace(2,10,33),-10.Alinspace(-5,0,15),10.Alinspace(-5,0,14), 

linspace(2,10,33)] 

[io— 4 , 10] 

VHI 


p: 1 

p: linspace(0.1,12,80) 


DSI 

Vo 

1 + 


p: 1 

p: linspace(0.1,12,80) 


GMLFI 

Vo [1 + a iv) q ) 

or. 10.Alinspace(-7,3,31) 

p: linspace(0.5,12,24) 

q: linspace(0.5,12,24) 

[-40,40] 

LPI 

(£)>£)' 

VO- 1 

p: [linspace(0.5,12,24), 10.Alinspace( 1.5,2,2)] 

q : [linspace(0.5,12,24), 10. Alinspace( 1.5,2,2)] 


CNDI 

Vo 

{1 + /3 cos[a=(</?-¥> 0 )]} 2 

vo- o 

ck : [linspace(0.1,1,30), 10.Alinspace(-3,-1.5,4), 10. Alinspace(0.5,2,4)] 

j3: [linspace(l,10,30),10.Alinspace(-2,-0.5,8),10.Alinspace(1.5,2,3), 

-10. A linspace(-2,1,7)] 



Tables: Parameter values expressed in Matlab code: linspace (a, b, n) means 
n equispaced values between a and b;10.Alinspace(a,b,n) means n log- 
equispaced values between 10“ and 10 6 . Parameter Vo and the reduced Planck 
mass Mpi always set to 1. The range of studied field values is ip £ [10 —4 ,40] 
unless otherwise indicated. 

The models, choice of parameter values and range of field values of Table [8] 
have been subjected to the numerical test described in subsection |4.1| for the several 
determinations of the spectral parameters n s ,a s , (3 S listed in Section[3] Let us sum 
up the conclusions of the most relevant cases: 

For the determination of spectral parameters PLANCK2013 combined with 
WP and BAO data using the A CDM + r + a s model (i.e. without running of the 
running) depicted in Table 1, the computation has concluded that all the models 
in Table [8] are disproved within a confidence limit 93%. The reason in all cases 
is that the parameter values that match the bound for r and the expected value of 
n s within 2 standard deviations can only provide values of the running a s that are 
too close to 0, far above the value a s = —0.021)to oio °f this determination, and 
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consistently in the same tail of the Gaussian distribution. 

For the determination of spectral parameters PLANCK2015 TT+lowP, without 
running of the running, of Table 5 the result changes drastically: the models 

1. Natural Inflation (NI) |[T4l 

2. Power Law Inflation (PLI) llT6ll 

3. Constant n s A Inflation (CNAI) l30l 

4. Constant n s B Inflation (CNBI) ll30l 

5. Open String Tachyonic Inflation (OSTI) f29l 

6. Generalised MSSM Inflation (GMSSMI) lt3TI 

7. Generalised Regularised Point Inflation (GRIPI) li32l 

8. Constant Spectrum Inflation (CSI) ll33l 

9. Constant n s C Inflation (CNCI) lf30l 

10. Dynamical Supersymmetric Inflation (DSI) 041 . 

are disproved with confidence at least 94.5%, and asking for the number N of e- 
folds of expansion that the model allows to be in the range [30, 80] only rules out 
completely a further model, namely Minimal Super-symmetric Standard Model 
Inflation (MSSMI) l35l . The rest of the models in the table admit some choice(s) 
of parameter and field values that simultaneously satisfies all filtering conditions. 

The inclusion of the running of the running in the numerical test results in 
another reversion of conclusions. For the determination of spectral parameters 
PLANCK2015 TT+lowP with running of the running of Table 6, the only models 
in Table [8] which are not disproved for any choice of parameter and field values 
with confidence at least 92.8% are: 

1. Kaller Moduli Inflation I (KMII) O- 

2. Kaller Moduli Inflation II (KMIII) ifTTl . 

3. Logamediate Inflation (LMI) lf36l . 

4. Twisted Inflation (TWI) ll37ll . 

5. Brane SUSY Breaking Inflation (BSUSYBI) f38l . 

6. Spontaneous Symmetry Breaking Inflation (SSBI) lf39l . 

7. Running-mass Inflation (RMI) fl40ll 
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8. Generalised Mixed Large Field Inflation (GMLFI) If34l . 

9. Constant n s D Inflation (CNDI) If33l . 

(albeit GMLFI is disproved with 92.1% confidence, and CNDI with 91.9% confid¬ 
ence). The most common reason for disproving the models is now that they only 
provide values of /3 S at more than 1.8 standard deviations of distance of the value 
0.029 of this determination, which we point out that is a quite large value compared 
with those of n s — 1 = —0.043 and a s = 0.011 in this determination. We will 
bring up this subject in the next section. 

Of the models that are not disproved for this PLANCK2015 TT+lowP with run¬ 
ning of the running determination of the spectral parameters, the e-fold test in our 
software has found only for the models KMIII.SBLSSBI,GMLFI,CNDI choices of 
parameter and field values passing the filters and further allowing the number of 
expansion e-folds N to lie in the range [30, 80]. 

5 Accuracy and reliability of the spectral parameter val¬ 
ues 

The results of the analytical and numerical fitting of the studied single field infla¬ 
tionary models to the successive sets of values provided by the Planck collaboration 
for the spectral parameters r, n s ,a s , f3 s may be summed up as: 

1. The spectral values provided by PLANCK2013 for r, n s ,a s , in models without 
running of the running, disfavour all of the studied models with a confidence 
limit in the ranges 93%-95.5%. The poor fit of the value of the running a s 
supplied by the models to the measured one, which has a larger than expected 
size, is the most common reason for this invalidation. 

2. The corrected values provided by PLANCK2015 for r,n s ,a s , in models 
without running of the running, may be fitted by most models in our study. 
The main reason for the change is the greatly diminished measured value of 
the running a s . 

3. But if one uses the values provided by PLANCK2015, for models including 
running of the running, for the parameters r, n s , a s , (3 S again most models 
are disfavoured with a confidence limit 92.8% or better. The most common 
reason is the poor fit these models provide for the running of the running j3 s , 
which in this measurement has a larger than expected size. 

The dramatic contrast in results leads to the question of the reliability of the 
successive determinations of the spectral parameters r, n s ,a s , f3 s . The authors have 
found that the methodology of model fitting used in Sl|,@ has a bias that often 
leads to the overestimation of the highest order parameter, which in turn triggers 
biases in the estimation of the lower order parameters. 
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The bias does not arise from the probabilistic algorithms, such as the Monte 
Carlo Markov Chain (MCMC) method, employed to determine the measured val¬ 
ues of the spectral parameters. It comes from the procedure to judge the fit of the 
model. 

It is another manifestation of the well known bias of the least square regression 
polynomial tit when the fitted polynomial p(x) has too small a degree: the highest 
order coefficient of p(x) is found to have a very large value, which is actually 
caused by the measured data growing faster than the degree of p(x) allows. The 
second highest order coefficient of p(x) typically has a bias of opposite sign to 
compensate for the overestimation of the highest degree term. Figure 1 illustrates 
this phenomenon with an elementary example. 

This problem is a particular case of the more general difficulty in asessing the 
value of a Taylor polynomial, that approximates a function in the neighbourhood 
of a Taylor expansion point, from values of the function at points that lie further 
and further from the expansion point. Conceding any weight at all to the values 
of the function at points far from the expansion point in the assessment of best fit 
carries a great risk of introducing a bias such as that described in Figure 1. 

The fitting of the model in every step of the MCMC algorithm of U,0, and 
in comparable algorithms, is probabilistic, but ultimately close to the classical re¬ 
gression fit: the spectral parameters on which the model depends are given values 
that maximize the likelyhood of the observations that have to be described by the 
model. The coefficients Ci of the observed power spectrum are assumed to follow 
a Gaussian distribution. This means that the spectral parameters determining the 
theoretical power spectrum V-jz(k) for scalar perturbations are selected in order to 
minimize the distance between the values Ci given by the model and the observed 
values Cj. 

This distance is expressed in standard deviations for the probabilistic fitting 
that maximizes likelyhood, and in natural units by the classical polynomial regres¬ 
sion fit that minimizes residue. But the variation of the standard deviation for each 
spectral parameter changes little as the parameter varies its value, so minimizing 
the distance in natural units or in standard deviations ends up assigning very sim¬ 
ilar values to the spectral parameters. Indeed, this is the reason why probabilistic 
methods such as MCMC are often preferred to regression analyses that are more 
comprehensive but far more computationally costly and end up reaching a similar 
result. 

A consequence of this equivalence of distances is that likelyhood-maximizing 
Bayesian, MCMC .. .methods, while computationally vastly more efficient than a 
regression fit, inherit the latter method’s bias, displayed in Figure 1: fitting a model 
for a function In V-jz(k) with parameters up to an insufficient degree will result in 
an overestimation of the highest order parameter, which in turn triggers a cascade 
of estimation errors in the lower order parameters. 

The order of the growth of the modelled function \nV-n{k) (the logarithm of 
the power spectrum) is currently unknown, but the pattern of systematically over- 
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estimating the highest order spectral coefficient in its Taylor series, and drastically 
revising it down as higher coefficients are incorporated to the model suggest that 
this bias is indeed happening: 

• The PLANCK2013 determinations of the value of the running a s in a A- 
CDM model without running of the running (table 5 of |[Q) attribute to a s 
expected values ranging from -0.0149 to -0.0094, 

• but the determinations by PLANCK2013 with the same methodology, adding 
a running of the running to the model put the expected value of a s in the 
range [0,0.006] (one order the magnitude smaller in absolute value than with 
only running; sign uncertain due to proximity to zero). 

• The value of the running of the running {3 S is estimated by PLANCK2013 
(table 5 of 0]|) to lie in the range [0.017,0.020], and by (19) in PLANCK2015 
to lie in the range [0.025,0.029]. These values arc one order of magnitude 
greater than those attributed to the running in these estimations. 

The values of the spectral parameters n s — 1, a s , (3 S will be reliably known only 
when higher order terms in the Taylor series of In 'P-r (/.:) are known, or a method 
without this bias is used to fit the models to the data. 


6 Conclusions 

Evidence mounts, both analytical and numerically, that single field slow rolling in- 
flaton models do not fit well the observations of the CMB radiation. Sophistications 
such as multiple fields, or a breakdown in the slow roll regime, or a completely dif¬ 
ferent paradigm such as the Matter Bounce Scenario ought to be contemplated. 

But the value of the spectral parameters n s — 1, a S: j3 s of the CMB radiation is 
not yet well established, so we are witnessing the twilight, rather than the death, of 
this family of models. 

In Section [3] of this work we have found analytically bounds for the values 
of the spectral parameters that most of the single field, slow roll inflation models 
support. In Section[4]we have described our MATLAB software package that finds, 
for each model given by its potential function V(p), all possible values of the 
spectral parameters n s — 1, a s , /3 S that the model supports, for a comprehensive list 
of values of the field ip and of further parameters on which the potential V(<p) may 
depend. There is a deliberate redundancy in the two approaches, and both have 
reached the same conclusions for every experimental determination of the values 
of the spectral parameters. The numerical testing software described in this work 
can be applied to any single field slow roll model and determination of the spectral 
parameters beyond those studied here. 

The conclusions reached by our analytic/numerical assessment of the likely- 
hood of single field, slow roll inflation models vary strongly with the different 
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determinations of the values of the spectral parameters, but there is a clear pattern 
in this variation: 

1. the most accurate determinations PLANCK2013 without running of the run¬ 
ning (such as PLANCK2013+WP+BAO) assign to the running a s a big 
value, that rules out most of the models because in them a s is a second 
order expression in the slow roll parameters and can only be much smaller 
in magnitude, 

2. the determinations PLANCK2015 give a much smaller value to the running 
a s , which can be fitted to most of the tested inflaton models, 

3. but if the value of the running of the running of PLANCK2015 is included 
in the likelyhood test, it turns out that its determination by PLANCK2015 
(both TT+lowP and TT,TE,EE+lowP) has such a big magnitude that again 
most of the tested inflaton models cannot furnish values within 1.8 standard 
deviations of the expected value. 

The reason for this pattern of contradictory conclusions seems to be a mathem¬ 
atical bias in the method that has been used for the computation of the values of 
the spectral parameters, which results in a systematic overestimation of the highest 
order one. The bias is a migration of, and very close to, a classical bias of regres¬ 
sion (i.e., mini miz ation of residue) fitting: if one tries to fit a polynomial of too low 
degree to values of a function that actually grows faster, no matter how flawless the 
procedure for finding the better fit is, it will result in a polynomial with an exag¬ 
gerate value for the magnitude of the leading term, which will result in a cascade 
of further errors for the lower coefficients as the table in Fig. 1 illustrates. 

The probabilistic (Bayesian, MCMC, ...) methods currently used to fit the 
value of spectral parameters following a Gaussian distribution in a model actually 
minimize the residue, expressed in standard deviations rather than in natural units, 
and reproduce this bias. 

The pattern of disproving/validating/disproving of the single field models by 
successive determinations of the spectral parameters is not completely symmet¬ 
rical, because the instances in which the models are validated are often extreme, 
narrow choices for the values of the parameters on which the model depends. 
Hence the authors’ suspicion that single field slow roll inflation models will ul¬ 
timately have to be discarded. 

Nevertheless, to rule out inflation models based on their fit to the measured 
spectral parameters r, n s ,a s , j3 s will not be possible until the value of these spectral 
parameters is reliably known, for which more terms in the Taylor series of In V-ji (k) 
or, even better, a determination procedure free of its current bias, will be required. 

This investigation has been supported in paid by MINECO (Spain), projects 
MTM2011-27739-C04-01, and MTM2012-38122-C03-01. 
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Figure 1: Data points on curve y = x A \ the regression polynomials ^ mx l are 
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The leading order coefficient in each regression polynomial is systematically over¬ 
estimated (from its true value 0) to try fitting the growth of the function, which is 
actually of a higher order. This bias cascades down to the lower order coefficients, 
starting typically with an underestimation of the second to highest order coefficient. 
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